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ABSTRACT 

We present evidence of Fe fluorescent emission in the Chandra hetgs spec- 
trum of the single G-type giant HR 9024 during a large flare. In analogy to solar 
X-ray observations, we interpret the observed Fe Ka line as being produced by 
illumination of the photosphere by ionizing coronal X-rays, in which case, for 
a given Fe photospheric abundance, its intensity depends on the height of the 
X-ray source. The hetgs observations, together with 3D Monte Carlo calcula- 
tions to model the fluorescence emission, are used to obtain a direct geometric 
constraint on the scale height of the flaring coronal plasma. We compute the Fe 
fluorescent emission induced by the emission of a single flaring coronal loop which 
well reproduces the observed X-ray temporal and spectral properties according 
to a detailed hydrodynamic modeling. The predicted Fe fluorescent emission is 
in good agreement with the observed value within observational uncertainties, 
pointing to a scale height < 0.3R*. Comparison of the HR 9024 flare with that 
recently observed on II Peg by Swift indicates the latter is consistent with exci- 
tation by X-ray photoionization. 

Subject headings: hydrodynamics — plasmas — stars: coronae — X-rays: stars 

1. Introduction 

Since the early 1970's, spatially resolved observations of the solar corona have revealed 
a high degree of structuring over all scales, from of order of the solar radius down to instru- 
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mental resolution limits (e.g., Vaiana et al. 1973al lbl). This spatial structuring is intimately 



linked to the characteristics of magnetic field generation and interaction, and to plasma 
heating mechanisms. On other stars, especially very active ones (with X-ray luminosity Lx 
up to 10 4 x Lx©), coronal structure and its relation with stellar parameters such as mass, 
rotation and evolutionary phase, remains very uncertain. 

Techniques used to d ate to investigate the morphology of stella r coronae comprise ro- 



Reale et al. 2004; Testa et al 



tational modulat i on (e.g.. iMarino et al.l 120031 ; iFlaccomio et al.l 120051 ) , flare m odeling (e.g., 



2007bl . hereafter Paper I), eclips e mapping (e.g. 



1990), spectroscopic density and radiation fie 
2004 ) , resonance scattering (e.g., 



velocity modulation (e.g. , iBrickhouse et al 



Testa et al 



d diagnostics (e.g.jTesta et al.ll2004a 



2001 



White et al. 



2004b; Matranga et al 



Chung et al 



2005 



2004 



Ness et al 



Testa et al. 



2007a), 



Hussain et al 



2005 



Huenemoerder et al. hood ), and simultaneous Doppler imaging and X-ray spectroscopy (IHussain et al 
2007h . 



The fluorescent iron line at ~6.4 keV presents a further, potentially powerful diag- 
nostic of stellar coronal geome try. Extensively used in study of active galactic nuclei and 
X-ray binaries (see e.g. reviews iGeoree fe FabiaJl99ll Ikevnolds fe NowalJ|20jd), and often 



seen on the Sun during flares (ICulhane et al 



1981 



Parmar et al. 



1984 



Tanaka et al. 1984 



Zarro et al.lll992l ; iPhillips et al.l 119941 ). the line is produced by electron cascade after one of 
the two K-shell (n—1) electrons of an iron atom (or ion) is ejected following photoelectric 
absorption of X-rays. In the solar and stellar case, for a given irradiating flare spectrum the 
line strength depends essentially on t hree parameters: the flare hei ght, heliocentric angle, 
and photospheric Fe abundance (e.g., lBailll979l ; iDrake et al.ll2007bl . hereafter D07b). If it 
can be observed in stars, Fe K fluorescence therefore presents a means of estimating flare 
and coronal scale height. 

Fe flu orescence has now been detected on a num ber of pre-main sequence stars with 
disks (e.g., iTsujimoto et al.ll2005l : iFavata et al.l 120051 ) , where it is most likely produced by 
coronal X-ray irradiation of the cold disk material. Here we present evidence for photospheric 
Fe fluorescent emission in the Chandra hetgs spectrum of the X-ray active single Gl giant 
HR 9024 (see Paper I and references therein for a description of the characteristics of the 
target). The observations and analysis are briefly described in §21 In §3] we present 3D 
Monte Carlo calculations of the Fe Ka fluorescence. The results are presented in §H where we 
derive an estimate for the coronal scale height and compare it with the prediction of the loop 



1 The Fe fluorescent feature actually consists of two components, at 6.391 and 6.404 keV (for Fei), which 
are unresolved by present instruments, therefore hereafter we will refer to the feature as to a single emission 
line. 
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hydrodynamic model that succeeds in describing the observed flaring spectrum (Paper I). To 
our knowledge, the HR 9024 Chandra spectrum represents only the second detection of Fe 
photospheric fluorescence. Very recently, the line was observed by the Swift X-ray Tel escope 
during a "superflare" on the RS CVn binary system II Peg by lOsten et al.l (120071 ). who 
ascribed the excitation mechanism to electron impact ionization of photospheric Fe instead 
of photoionization. We discuss our result and the II Peg observations in £j5j and show that 
the latter is also consistent with, and more plausibly associated with, photoionization than 
electron impact. 



2. Observations and Data Analysis 



Chandra High Energy Transmission Grating (ICanizares et al.l 120051 ) observations of 
HR 9024 (ObsID 1892) were discussed in detail by Paper I that analyzed the X-ray spectra 
and variability of HR 9024 and constructed a flaring loop hydrodynamic model that success- 
fully reproduce d the observations. We analyse the same data here, using the PLNTofALrf] 
IDlJf software ( Kashvap fe Drake 2000 ), and t he Interactive Spectral Interpretation System 
(ISIS0) version 1.4.2 iHouck fe Denicolalboooh . 



The heg flare spectrum (first 40 ks of the observation) in the wavelength region with the 
Fe Ka emission is illustrated in Figured] (top panel), together with a comparison between the 
evolution of the fluorescent and adjacent continuum emission {bottom panel). A Kolmogorov- 
Smirnov test applied to the lightcurves of Fe Ka and the nearby continuum indicates that 
they are significantly different during the flare (P < 10~ 10 ), but P ~ 0.94 during the rest of 
the observation. The Fe Ka line flux was measured with ISIS by convolving a Gaussian profile 
with the instrumental response, using the continuum predicted by an isothermal spectral 
model described in §|3] (and derived in Paper I). We obtained an Fe Ka line flux of 1.7 x 
10~ 5 ph cm -2 s _1 , with a corresponding 90% confidence range 0.8 — 2.8 x 10 -5 ph cm -2 s _1 . 



3. Fluorescence Modeling 



The basic physics of Fe fluorescent line production on the Sun has been investigated 
through detailed Monte Carlo calculations of the interaction between a hard X-ray continuum 



2 http: / /hea- www.harvard.edu/PINTofALE 
interactive Data Language, Research Systems Inc. 
4 ISIS is available at http://space.mit.edu/cxc/isis/ 
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(E > 7.11 keV) and dense, neutral photospheric material by lBailll979l . considering a thermal 
continuum and taking into account the curvature of the solar photosphere. The problem has 
been further discussed in detail in the stellar context by D07b. 



We used a modified version of th e 3D radiative transfer code MOCASSIN ( lErcolano et al. 



20031 . 120071 ; iDrake fe Ercolanoll2007l ). to carry out Monte Carlo calculations to follow the Fe 
Ka photons resulting from X-ray absorption and to determine the emerging Fe Ka emission 
as a function of the flare geo metry par ameters illustrated in Figure [21 Here, our method 
differs from th e calculatio ns of lBail (119791 ) in that we can use any given input fluorescing spec- 
trum, whereas iBail (119791 ) adopted an analytical expression for the thermal bremsstrahlung 
continuum and computed a grid for T in the range 0.5-5 keV. Detailed descriptions of the 
methods used for the fluorescence calculations are presented by D07b. 

We note that the fluorescence efficiency, £ , defined as the ratio of the Fe Ka flux to 
one-half the integral of the X-ray flux above 7.11 keV (£ = I(FeKa)/[I c (E > 7.11 keV)/2]) 
is only very weakly dependent on the photospheric metallicity, [M/H] (D07b). For HR 9024 
no accurate photospheric abundances were available, and we therefore obtained optical spec- 
tra with the high resolution (R = 86000) spectrograph SARG at the Telescopio Nazionale 
Galileo, and derived a near-solar photospheric Fe abundance [Fe/H] = —0.2 ± 0.3 (L. Affer 
et al. in prep.). Based on the calculations of D07b, we find an uncertainty in S of ±25% for a 
metallicity uncertainty range [M/H]= —0.5 — hO.l, including allowance for an o verabundance 
of q- elements with respect to Fe at the low metallicity end of 0.2-0.3 dex (e.g. 
2003h . 



Bensby et al. 



In order to compute £ we used a fluorescing spectral model with logT[K] = 7.8, and 
logEMfcm" 3 ] = 54.9, based on the values found from the analysis of the continuum emission 
in the early phases of the flare (see Fig. 6 of Paper I). This was preferred over a spectrum 
computed using the emission measure distribution, EM(T), derived for the flare as a whole 
(rise + peak + decay; see Figs. 3 and 6 of Paper I) because most of the Fe Ka emission is 
detected in the rise, peak, and early decay phases: the use of the whole flare EM(T) would 
slightly underestimate the hard X-ray flux. Values of £ were calculated assuming a point 
source for the flare and for a range of flare heights, h. The results obtained for a range of 
flare heights and line-of-sight viewing angles are shown in Figure [3J 

As expected, the model fluorescent efficiency declines strongly for larger values of h, the 
height of the flaring source. This occurs because of the 1/h 2 dilution of the irradiating surface 
flux, and because photon incidence angles become larger, causing the average depth where 
the photoionization takes place to increase. For the line-of-sight viewing angle <fi < 60°, the 
fluorescence efficiency depends primarily upon h but not on <ft; for cf) > 60°, the efficiency 
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rapidly decreases^. 

In Paper I we modeled the observed flare using a hydrodynamic treatment of the evo- 
lution of the flaring plasma within the confines of a single long loop. The hydrodynamic 
model provides us with the physical conditions (in particular T and n e ) of each plasma vol- 
ume element along the loop, allowing us to compute the expected Fe Ka emission for a direct 
comparison with the measured fluorescence efficiency In order to calculate the £ induced by 
the emission of the flaring loop as resulting from the hydrodynamic modeling we proceeded 
as follows: (1) we synthesized hard X-ray spectra (7-50 keV) for isothermal plasmas for a grid 
of temperatures covering a large T range (logT = 5, 8.5); (2) we divided the half-loop in four 
different segments and computed hard X-ray spectrum of each segment by adding the spectra 
of all volume elements (each characterized by a temperature and emission measure) within 
the segment, integrating over the 40 ks of the flare; (3) we computed the expected fluores- 
cence emission induced by each segment, depending on its input spectrum and corresponding 
height, and computed the overall Fe Ka emission, and therefore the predicted £. Splitting 
the loop in segments linearly spaced in h (i.e., hi = {0.125, 0.375, 0.625, 0.875} x h), assuming 
that the loop is perpendicular to the stellar surface, we obtain £(HD mod) = 0.016. We also 
computed the efficiency using a logarithmic spacing (hi = {0.0065,0.028,0.118,0.504} x h), 
and we find £(HD mod) = 0.018 only slightly different from the former value. 



4. Results 

To compare the observed Fe Ka flux with model values, we used the same spectral model 
based on the observed continuum described above. The one-half integrated X-ray flux above 
7.11 keV (up to 50 keV) as computed from this model is 7.5 x 10~ 4 ph cm -2 s _1 , which, 
when combined with the observed Fe Ka flux, yields a value for the observed fluorescence 
efficiency £ = 0.022 [0.010 — 0.033]. We note that at 50 keV the flux is already several orders 
of magnitude lower than at 7.11 keV, therefore the efficiency does not depend significantly 
on this upper bound value used for the energy range. In Figure E] the measured value, with 
the 90% confidence range, is represented as a shaded area superimposed to the theoretical 
curves derived in §|3j 

The comparison between observed and mod el values of £ imp lies for the scale height 



an upper limit of ~ 0.3-R*, where R+= 13.6i? (jSingh et al.l 119961 ). This result compares 
well with the hydrodynamic modeling that finds a loop semi-length L ~ 0.5/2* (Paper I), 
corresponding to a loop apex at h = L x 2/n = R*/ir ~ 0.32/2* (assuming that the loop is 



The results are qualitatively similar also for the case of the semi- infinite slab in 



George fe Fabian 



(|1991I ). 
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perpendicular to the stellar surface, otherwise h < 0.32i?*). The efficiency expected on the 
basis of the hydrodynamic model itself corresponds to an "effective scale height" h ~ 0.1R+— 
somewhat lower that the maximum height of the loop model. This can be understood on the 
basis of the extreme characteristics of the flaring loop in HR 9024, which is much hotter than 
typical solar flares: for solar flares, typically reaching temperatures > 10 7 K, the hard X-ray 
photons above the 7.11 keV threshold will be emitted almost uniquely by the hotter loop 
apex; instead, in the flaring loop on HR 9024, with peak temperature of the order of 10 8 K, 
the lower coronal plasma will still be effective in producing hard X-ray photons because of 
its elevated temperature (of several 10 7 K) and density. 



Discussion and conclusions 



We have presented the first detailed analysis of observed photospheric fluorescence to 
probe the geometric properties of stellar coronal emission. For the flare on HR 9024 we find 
the observed Fe Ka intensity completely consistent with production through photospheric 
irradiation by flare X-rays. Our estimate for the flare scale height of h < 0.3-R* provides a 
cross-check for the results of the hydrodynamic modeling of the flare observed on HR 9024 
that can be satisfactorily reproduced by a single flaring loop of semi-length L = 0.5R*. The 
agreement between the observed 8 and the value computed directly from this model provides 
further confidence both in the hydrodynamic modeling approach and X-ray fluorescence as 
the excitation mechanism for the Fe Ka line. We note that our model based on Chandra 
data might underestimate the high energy continuum flux because we cannot detect any 
additional (either thermal or non-thermal) emission possibly present at energies outside the 
Chandra band. However, if additional high energy emission were present, the higher flux 
above 7.11 keV would make the efficiency lower bringing it in even better agreement with 
the predictions of the model. 

As noted in jjH Fe Ka emission has been extensively observed in solar X-ray spectra 
during flares 



Phillips et al. 



Culhane et al. 1981; Parmar et al. 1984: Tanakaetal. 1984; Zarro et al. 



1992 



1994 ). Two main production mechanisms were initially explored to explain the 



line: electron impact from a non-thermal electron population, and X-ray photoionization. 
The observed Fe Ka evolution with respect to the thermal and non-thermal emission, and the 
observed ce nter-to-limb variation s seen on the Sun, s t rongly support the X-ray fluorescence 
mechanism (ICulhane et al.l 11981c iParmar et al.l 1 1 9 84j ; iTanaka et al.l Il984l ) , although in rare 
cases there is indication that e lec tron impact might contribu te during the early impulsive 



phase (e.g., lEmslie et al.lll986l ). iBallantyne fe Fabian! (120031 ) have also shown that Fe K 



production in accretion disks by non-thermal electron bombardment is extremely unlikely, 



- 7- 



and requires 2-4 orders of magnitude greater energy dissipation in the electron beam than 
is required for an X-ray photoionization source. In this context, the rece nt observations of 



an ex tremely large Fe Ka equivalent width in the PMS star V1486 Ori (jCzesla & Schmitt 



20071 ). a nd variability in Fe Ka emission uncorrelated with the observed X-ray continuum on 
Elias 29 (jGiardino et al.ll2007l ). seem to challenge the photoionization excitation mechanism. 
However, D07b note that fluorescence from PMS disks excited by X-rays originating from 
the unseen stellar hemisphere will also be observed: in < 50% of cases, Fe Ka can then 
appear anomalous when compared only to the observed thermal continuum. 



To date, the Swift observation of II Peg (lOsten et al.ll2007l ) is the only other case where 
Fe Ka emission has been obse r ved fo r a star (other than the Sun) lacking substantial circum- 
stellar material. lOsten et al.l (120071 ) dismissed X-ray fluorescence as the excitation mecha- 
nism based on arguments for fluorescence in an optically-thin medium, and instead attributed 
the Fe Ka to electron impact with non-thermal el ectrons. As noted in earlier discussions 



(e.g. iParmar et al.lll984l ; iBallantyne fe Fabianll2003l ). one major problem for an electron im- 
pact excitation mechanism is the inefficiency of this process. Extremely large energies must 
be dissipated in accelerated electron beams to produce continu um and fluorescence emission 
that might be observed on a star. Indeed, lOsten et al.l (120071 ) showed that the hard X-ray 
continuum observed in the II Peg flare would require an energy of 3 x 10 40 ergs in non-thermal 
electrons if interpreted in terms of thick-target bremsstrahlung. Over the time of the dura- 
tion of the flare, this corresponds to an energy dissipation rate in accelerated electrons of 
more than 100 times the stellar bolometric luminosity. The hard X-rays are also interpreted 
as lasting for timescales comparable with the flare soft X-ray emission, at variance with the 
Sun where this non-thermal component, when present, is generally impulsive. Regardless of 
the true nature of the hard X-ray flux, the equivalent width measured from the II Peg Swift 
spectra ranges between 18 and 60 eV, corresponding to efficiency values £ between 1% and 
2%. These values are well within the range found from the theoretical calculations presented 
here, and comparable with the measured value for HR 9024. Since the non-thermal hard 
X-ray component in the II Peg flare is not well-constrained by the Swift data and could also 
be explained by thermal emission, the observed Fe Ka line seems more easily explained as 
arising from photoionization by flare X-rays (see also the discussion of D07b). 
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Fig. 1. — Top: HEG spectrum of HR 9024 during the flare (first 40 ks of the observation) 
in the wavelength range containing the Fexxv coronal emission (~ 1.85A) and the Fe Ka 
emission (~ 1.94A), with the best fit model. Bottom: Comparison of lightcurve of the Ka 
emission (black) with the adjacent continuum background emission (red); the continuum is 
extracted on a wavelength region (1.912-1.932A+ 1.948-1.968A) that is 2.5 times larger than 
the wavelength extraction region used for the Fe Ka, and then appropriately scaled. We also 
plot the scaled lightcurve of dispersed photons integrated in the HEG band (green curve). A 
temporal binsize of 5 ks is used for all lightcurves. We plot counts (scaled by 0.003 in the 
case of the integrated HEG band), therefore the associated uncertainties correspond to the 
Poisson errors. 
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Fig. 2. — Geometry for fluorescent excitation of photospheric Fe Kq by coronal X-ray emitted 
by a source located at height h above the photosphere. <p is the angle of inclination with 
respect to the line-of-sight. 
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Fig. 3. — Theoretical curves of fluorescence efficiency (black curves), obtained using the 
MOCASSIN code, as a function of the angle of inclination with respect to the line-of-sight, 
4>, for different heights of the hard X-ray source: h = 0.001-R*, 0.1R+, 0.5-R*. We superimpose 
the measured value (blue dotted line) with the corresponding 90% confidence range (shaded 
light blue area), and the prediction of the hydrodynamic loop model (red shaded area). 



